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Abstract — The primary objective of this paper is to find Pareto optimal solutions to 
multiobjective linear programming problem with fuzzy coefficients under imprecise 
environment. There exist several methods of Pareto optimization in literature. Recently 
Jimenez, and Bilbao showed that fuzzy efficient solution may not be Pareto optimal when 
one of the fuzzy goals is fully achieved. Next Wu, Liu, and Lur redefined membership 
function of fuzzy set theory. But we have observed that redefined membership function was 
not used during computation according to that proposed definition. And proper usage of 
membership or even redefined membership functions in fuzzy environment may suggest 
problem as infeasible. Moreover in decision making under imprecise environment, 
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membership function is constructed so that higher membership grade gives better result 
for objective. And it is sometimes difficult or even impossible to construct fuzzy 
membership function or redefined membership function of imprecise objectives. And usage 
of membership function or redefined membership function in fuzzy environment gives at 
most one Pareto optimal solution to the problem. So in this paper, functions, namely T- 
characteristic functions, are used in place of usual membership functions or redefined 
membership functions in fuzzy environment. It results in more than one Pareto optimal 
solutions, and ultimately this feature helps decision maker. We have proposed a 
mathematical procedure to further illustrate usage of T-characteristic function in finding 
Pareto optimal solution to multiobjective linear programming problem with fuzzy 
coefficients under imprecise environment. Numerical examples illustrate our proposed 
procedure. 

Index Terms — Fuzzy set theory, Linear programming, Optimization, Pareto 
optimization, T-characteristic function. 

I. INTRODUCTION 

We consider the following multiobjective linear programming problem (MOLPP) with k 
objective functions Z t (x) = c t x , i = l. . . k as 

Min z(x) = (z 1 (x),z 2 (x),...,z k (x)) T =(c 1 x,c 2 x,...,c k x) T 

(A) 

subject toxsX 
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HereX ={xe R n :Ax<b,x>0},b = (b 1 ,b 2 ,...,b m )e R m ], and A is a mxnmatrix. R denotes the set 
of reals. By assuming that the decision maker (DM) has imprecise aspiration levels for each of 
the objective functions z, ■(*) = c t x , i = 1. . .k of model (A), researchers has proposed several 
methods in literature for characterizing Pareto optimal solutions to MOLPP (A). In one such 
approach, we use fuzzy set theory [7]. Here we convert MOLPP (A) into single objective 
optimization problem to obtain M-Pareto optimal solution, and then Pareto optimal solution by 
using the fuzziness of the DM's aspiration with respect to the goals of imprecise objective 
functions (constraints till they are symmetric). 

One of the popular, well established methods to solve fuzzy MOLPP is two phase method. 
Here max-min operator is usually applied due to its easy computation but it does not guarantee to 
yield an M-Pareto optimal solution if more than one optimal solution exists. 

Recently Jimenez, and Bilbao (2009) [9] observed that a fuzzy-efficient solution may not be a 
Pareto optimal solution in the case that one of fuzzy goals is fully achieved. And they proposed a 
procedure which extended the two-phase approach of Guu, Wu (1997, 1999) [14] and approach 
of Dubois, Fortemps (1999) [4] to attain a Pareto-optimal solution for model (A). 

But according to Yan-Kuen Wu et al (2015) the proposed approach by Jimenez and Bilbao 
(2009) [9] cannot guarantee to be a general procedure to attain Pareto-optimal solution to model 
(A). In their proposed approach, Wu et al have modified membership function of fuzzy set 
theory and introduced redefined membership function. But that definition of redefined 
membership function /^-(c-*) and its usage in mathematical models as well as given numerical 
examples in their paper are not analogous to one another. From mathematical model or from 
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given numerical examples in that paper, it is clear that only one part of redefined membership 
function (c-x) where value of objective function does not exceed sum of goal, and tolerance 
(for minimization type of objective functions) is used. Yet these additional constraints were not 
added to set of existing constraints. Thus in problem formulation redefined membership function 
Mi (.CfX) is strictly monotone over entire domain where as it is not so in definition. 

Actually Wu et al (2015) [17] removed upper bound of membership function only but there 
still lies a lower bound of it at zero. In those examples given by Wu et al, those additional 
constraints in definition of redefined membership function that objective value cannot exceed 
sum of goal, and tolerance (for minimization type of objective functions) have made no impact. 
But in many cases, if the lower bound of membership function remains intact, the problem may 
become infeasible. Point to be noted. 

Deep analysis highlights the cause to lie in the definition of membership function of fuzzy set 
itself. Membership function of fuzzy set provides satisficing result when extreme ends of 
imprecise information can be quantified within the boundary of zero, and one. But it may not be 
always logically correct to quantify or measure every type of imprecise information with in some 
bounded subset of real line. 

It is noted that an element of a fuzzy set can lie partly or never lie or must lie in that set. The 
idea of lying, lying partly or not lying is well quantified / measured by membership function of 
fuzzy set theory. But it does not suit in case of must lying (membership value always one), and 
never lying (membership value always zero). So usage of another function (in place of 
membership function or redefined membership function of fuzzy set) that best describes the 
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underlying issues well is now necessary. 

Moreover in the models that closely describe as well as represent real life situations, the 
objective functions, and the constraints in MOLPP may involve parameters or coefficients many 
of whose possible values are usually assigned by the DM. In the conventional approach, those 
parameters are given fixed crisp values in an experimental and/or subjective manner through the 
understanding of DM. But in real life situations, however, it is natural to consider that the 
possible values of these parameters are often only ambiguously known to the DM. In this case, it 
becomes more appropriate to use fuzzy numbers as values of parameters. The resulting 
optimization problem involving fuzzy coefficients may be viewed as more realistic, and 
preferred version than the conventional one. Historically Tanaka, and Asai formulated the 
MOLPP with fuzzy parameters [7]. Sakawa (1986) used comprehensive, interactive technique to 
find fuzzy satisficing solution of MOLPP with fuzzy parameters [10]. Various other researchers 
have also worked on it in recent years. 

We arrange this paper as follows: We first discuss the problem as well as existing methods to 
find solutions of MOLPP with fuzzy coefficients. The problems of using classical membership 
function or redefined membership function lead us to T-characteristic function [1], and 
consequently parametric T-characteristic function. Next T- a -Pareto optimal solution is defined 
analogous to existing M- a -Pareto optimal solution. Then we convert MOLPP with fuzzy 
coefficients under imprecise environment to a single objective linear programming problem by 
using proposed parametric T-characteristic functions, and discuss related lemmas. We integrate 
all those and develop an algorithm to find a -Pareto optimal solution to a -MOLPP for fixed a 
by using suitable parametric T-characteristic functions. First numerical example shows 
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advantages of using parametric T-characteristic function in place of membership or redefined 
membership function to solve a -MOLPP with fuzzy coefficients whereas second example is 
based on Eatman et al (1979), and shows how we generate different a -Pareto optimal solutions 
using different set of parametric T-characteristic functions under imprecise environment. 
Conclusions are drawn at last. 

II. Overview of Multiobjective Linear Programming Problem with Fuzzy Coefficients 
Let us consider the following MOLPP with fuzzy parameters having k objective functions 
Z i (x) = C i x, i = 1. . . k as 

Min z(x) = (C^,C,x C , x) T 

k (1) 

Subject to x e X(A,B). 
HereZ(A,B) = {xe R" : A J x<B j ,x>0,j = L..m},and C ( =(C,. 1 ,C,. 2 ...C,.„) is an n-dimensional row 
vector of fuzzy coefficients involved in the ith objective function c ( x, A. = (A n ,A j2 ... A jn ) is an 
n dimensional row vectors of fuzzy coefficients involved in the jth constraint AjX < B j , andi^ is 
a fuzzy coefficient of right hand side constant in the jth constraint A.x < B ] . These coefficients 
reflect the DM's ambiguous understanding of the nature of the coefficients involved in the 
process of problem formulation, and hence all these coefficients 
C ;1 ,C i2 ... C in ,i = 1... k;A n ,A J2 ... A jn , j = 1, 2 ... m and B } , j = 1... m in the MOLPP with fuzzy 
coefficients are fuzzy numbers with membership functions n c (C n ),{i c (C i2 )...Mc. (C in ),i = l...k; 
jU A i (a jl ), J u A 2 (a j2 )...jU A (a jn ) j = 1,2... m; /J B ib j )J=l..m. For the sake of simplicity alone, the 
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following notations may be used c = (c i ...c k ),a = (a 1 ...a m ),b = (b i ...b m ),and 
C = {C l ...C k ),A = {\...A m ),B = {B l ...B m ). 

Here fuzzy numbers may be involved in both objective functions, and constraints in the 
MOLPP with fuzzy coefficients. And so the notion of Pareto optimality cannot be applied 
directly and it requires extension. 

Again or -level set(0<a<l)of the fuzzy numbers A jr ,B p C ir is defined as crisp set(A, B, C) K 

for which the degree of membership functions exceed the level a (0 < a < 1) i.e. 

(A, B, C) a = {(a,b,c) : fi^ (a jr ) > a,jU Bj ty) > a,ju q (c ir ) > a; 
i = \...k; j = l...m; r = \...n). 

Suppose DM decides that degree of all membership functions of fuzzy numbers involved in the 
MOLPP with fuzzy parameters should be greater than or equal to some level a . So for such a 
given level a (0 < a < 1) , we select the coefficient vector (a, £>,c)e (A,B,C) a . Therefore the 
MOLPP with fuzzy parameters reduces to non-fuzzy MOLPP that depends on the coefficient 
vector (a, b,c)e (A,B,C) a in the form 

Min z(x) = (CjX,c 2 x... c k x) T 

subject to (1A) 
xg X(a,b)={xe R n : a .x < bj, x > 0, j=l...m} 

In fact infinite numbers of such MOLPP with fuzzy parameters exist for different values of 

coefficient vector (a, b,c)e (A,B,C) a , and (a, b,c) are chosen arbitrarily in the sense that degrees 

of all the membership functions for the fuzzy numbers exceed the level a (0 < a < 1) . Actually the 

DM desires to choose (a,b,c)e (A,B,C) a so as to minimize the objective functions under 
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constraints for given a (0 < a < 1) . Hence it is natural to formulate non-fuzzy MOLPP from 
MOLPP with fuzzy coefficients, and a certain level a (0 < a < 1) as follows: 
Min z(x) = (c l x,c 2 x,...,c k x) T 

subject to (2) 
xe X(a,b)={xe 7T:a.x<b.,x>0,j=l...m},(a,b,c)e (A,B,C) a . 

As per existing literature, it is a -MOLPP. Thea-Pareto optimal solution, and then a -level 

optimal parameters can be obtained by regarding (x, a, b, c) as the decision variables in the a- 

MOLPP. 

III. DEFINITIONS 

In the case of decision making under fuzzy environment [13], one interesting, and useful 
property of membership function of fuzzy set is that higher value of membership function always 
gives better result for objective function. And there may exist not only membership function but 
also many such functions with this amazing characteristic. Hence one new function viz. T- 
characteristic function that handles underlying issues as well as preserve this amazing 
characteristic of membership function well, is defined. 

T -characteristic function: Let S denotes the universal set. Then T-characteristic function of 
fuzzy subset A of S is denoted by T A , and is defined as T- A : S — > R . It assigns a real number 
T- A {x) to each element xe S , where T~ A {x) represents the grade of membership of xe S in A. Here 
R denotes the set of real numbers. 

Parametric T-characteristic function: Let S denotes the universal set. Then parametric T- 
characteristic function of fuzzy subset A of S with parameter a(0 < a < 1) is denoted byT A (x,a) , 
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and is defined asT A :S^R. For each a(0 < a < 1) , it assigns a real number T A (x,a)to each 
element xe S , where T A (x,a) represents the grade of membership of xe S in A. 

Next, on the basis of a- level sets of fuzzy numbers, the concept of <2- Pareto optimal 
solution to the <2- MOLPP (2) as a natural extension of the Pareto optimality concept for the 
MOLPP (A) is given as follows: 

a -Pareto optimal solution: A point x*e X(a*,b*) is said to be aa- Pareto optimal solution 
to a a -MOLPP if and only if there does not exist another xe X(a,b) , and {a,b,c)e (A,B,C) a 
such that c { x< qx*,i = l,...,£with strict inequality holding for at least one 1 . The corresponding 
values of the parameters (a*,b*,c*) are called a -level optimal parameters. 

Based on T-characteristic function, we now introduce the concept of T- a -Pareto optimal 
solution which is defined in terms of T-characteristic function instead of objective function. 

T- a -Pareto optimal solution: A point x*e X(a*,b*) is said to be a T- a -Pareto optimal 
solution to an a -MOLPP if and only if there does not exist another xe X (a, b) , and 
(a,b,c)e(A,B,C) a such that Tic^^^Tic^a) for all i=l...k, i # j, and 
T. (CjX*, a) < Tj (CjX, a) for at least one j for some fixed a (0 < a < 1) . The corresponding values of 
the parameters (a*, b*, c*) are called a -level optimal parameters for same a (0 < a < 1) . 

E.g. in MOLPP with fuzzy coefficients under imprecise environment, in case of maximization 
type of objective function zi; one of parametric T-characteristic functions may be zi itself. And 
in case of minimization type of objective function z 2 ; one of parametric T-characteristic 
functions may be (-) z 2 . They satisfy the same strictly monotonic property that higher value of 
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membership function always gives better result for objective function. And these may be used to 
convert a MOLPP with fuzzy coefficients or a -MOLPP into single objective optimization 
problem, and obtain T- a -Pareto optimal solution and hence a -Pareto optimal solution for given 
a (0<or<l). 

We note that the concept of T- a -Pareto optimal solution defined in terms of T-characteristic 
function of fuzzy set is a natural extension of M- a -Pareto optimal solution defined in terms of 
membership function of fuzzy set which is another natural extension of a -Pareto optimal 
solution defined in terms of objective functions. 

IV. Obtaining T- a -Pareto Optimal Solution and a -Pareto Optimal Solution to 
MOLPP with Fuzzy Coefficients 

Within the scope of the multiobjective decision making theory, the Pareto-optimality of 
solution is a necessary condition in order to guarantee the rationality of a decision. In a 
minimization problem, a fuzzy goal of DM may be to achieve "substantially less" than some 
value £. This type of statement can be quantified by eliciting corresponding suitable parametric 
T-characteristic functions. In some cases such description about goal may be completely 
unavailable. Then other techniques may also be employed to construct parametric T- 
characteristic functions of objective functions under fuzzy environment. 

Using the method as discussed by Bellman, and Zadeh, for fixeda,(0<a<l), (2) is equivalent 

to 
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max (T l (c l x),T 2 (c 2 x)...T k (c k x)) 

subject to (2A) 
xe X(a,b)= {xe R n : a jx < b ■ , x > 0, j = 1 ... m},(a,b,c) e (A,B,C) a . 

And this model (2A) is equivalent to 

Min {Cq (c^))" 1 , (T 2 (c 2 x))~ l ...{T k (c^x)))" 1 } T 

subject to (2B) 
xe X(a,b)={xe R n :&jx < bj,x > 0,j = l...m},(a,b,c) e (A,B,C) flr . 

Using the min-max method by Bowman (1976) [16], problem (2B) may be rewritten as 

Min max{(7} (c-x)) _ 1 }, i = 1, 2...k 

subject to (2C) 
xs X(a,b)={x g R n : &jX < bj, x > 0, j = l...m}, (a,b,c)e (A,B,C) a . 

Let v =max{(7:(c,.x))" 1 },/ = l,2...k . Then model (2C) reduces to 

min v 

subject to 

(7J (c-x)) _1 < v, ? = 1, 2 . . . k, ( 2D ) 
xe X(a,b)={xe /?":a ; -x <b y -,x >0,j = l...m}, (a,b,c) e (A,B,Q a . 

If we put v = J/j\then for same fixed or (0 < a < 1), it further reduces to 

max A 
subject to 

T i (c i x)>X,i = \,2...k, (3) 
xe X(a,b)={xe 7?":a ; -x < b ; -,x > 0,j = l...m}, (a,b,c)e (A,B,C) a . 

The relationships between the T-«r-Pareto optimal solution x* of the problem (3), and the a- 
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Pareto optimal solution to the generalized MOLPP with fuzzy coefficients for fixed a i.e. 0C- 
MOLPP (2) can be characterized by following lemmas. 

Lemma 1: If x*e X(a^,b R ) is unique T-CC -Pareto optimal solution to the max-min problem 
(3), then x* is a -Pareto optimal solution to the generalized a -MOLPP (2) for fixed a (0 < a< 1) . 

Proof: Let x*e X(a^,b R ) is unique T- a -Pareto optimal solution to problem (3) for fixed 
a (0 < a < 1) . If possible let x* is not a a -Pareto optimal solution to (2). Then there exist at least 
one yeX(a.a,b*), and (a,b,c)e (A3,Q a such that c t y < c,x*for all i, i^j, and c ^ < CjX* for at least 
one j. Since each7^(c,x), i = l,.„,Hs strictly monotonic decreasing function for each objective 
functions z i = c t x, i= 1, 2...k of minimizing type, we get by using this relation T t (c ( jc*) < T t (c,.y) 
V i =1, 2...k , i * j and Tj(CjX*) < Tfay) . This gives ye X{a L a ,b R a ) , and(a,b,c)e (A,B,Q a as 
an optimal solution to (3). It contradicts that x*e X(a^,b R ) is unique optimal solution to (3). 
Therefore our assumption is wrong. Hence x*e X(a^,b^ ) is a-Pareto optimal solution to a- 
MOLPP (2) for same a (0 < a < 1) ■ 

Lemma 2: If x*e X(a^,b R ) is a -Pareto optimal solution to the generalized a -MOLPP (2), 
then x* is a T- a -Pareto optimal solution to the problem (3) for some 7), i = l...jfcfor fixed 
a(0<a<l). 

Proof: Let x* e X (a% , b§ ) is a -Pareto optimal solution to MOLPP (2) for some a (0 < a < 1) . If 
possible let x*e X(a„,b*) is not a T- a -Pareto optimum solution to the problem (3) for any T l , i 
= 1, 2...k forar(0<a<l). Then for T n i = 1... k, there exists jce X(a^,b R ) , and (a,b,c) e (A,B,C) a , 
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we haver.(c.x) > r.(c.x*) = /l(x*),for somej, and^c,.*) > Tfcx*) = A(x*\ i = l...k,i * j . But each 
T; is strictly decreasing function for i = l...k. Hence we get CjX<CjX* for some j, and 
c t x< c t x*,i = 1,2,.. .k,i ^ j. It contradicts that x* is <2f- Pareto optimal solution to the <2- MOLPP 
(2) for#(0<a<l). Therefore our assumption is wrong. Hence x* is an optimal solution to the 
problem (3) for T i , i =1 . . . k for same a (0 < a < 1) ■ 

From lemma 1 , lemma 2, it is clear if the uniqueness of the optimal solution x*e X(a^,b^) , 
and(a,b,c)e (A,B,C) Q .for the max-min problem (3) is not guaranteed, it is necessary to perform a 
-Pareto optimality test of x*. And the «r-Pareto optimality test for x* can be performed by 
solving the following linear programming problem with the decision variables 
x-(x l ,x 2 ...x n ) T and £ = (£ l ,£ 2 ...£ k ) T . 

max^ e i 

Subject to (4) 

r,.(c,.x) - £, > r,.(c,x*),/ = i... k, 

x e X(a£,b*),(a,b,c)e (A,B,C) a ,e i >0,i = l...k 
Lemma 3: Let x and £ be optimal solution to (4) for some a (0 < a < 1) . Then 

1) If all £ i = 0 , then x* is a -Pareto optimal solution of the a -MOLPP. 

2) If at least one £ i > 0 , then x* is not a -Pareto optimal solution of the a -MOLPP. 
Instead of x*, x is a -Pareto optimal solution to the a -MOLPP. 

Proof of similar lemma is recorded in literature. It may be noted that several authors have 
proposed different other procedures to test a -Pareto optimality and can be used in place of (4). 
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V. PROCEDURE TO SOLVE MULTIOBJECTIVE LINEAR PROGRAMMING PROBLEMS WITH FUZZY 
PARAMETERS UNDER IMPRECISE ENVIRONMENT 

The above ideas can be further integrated into a general framework and develop an algorithm 
to obtain a -Pareto-optimal solution to a -MOLPP with fuzzy coefficients for given a (0 < a < 1), 
which is also T- a -Pareto optimal solution. The steps of the proposed algorithm are as follows: 

1) Elicit suitable parametric T-characteristic function ^(qx) for each imprecise objective 
function in such a way that higher value of parametric T-characteristic function gives better 
value of objective function. Well-posed parametric T-characteristic function preserves this 
characteristic and gives finite value. 

2) DM selects initial value of the parametera(0<a<l). Otherwise, some suitable value of a 
can be taken. 

3) Construct the max-min problem as (3) for chosen a(0 < a < 1) . 

4) For this fixed a(0 < a < 1) , solve that max-min problem, and obtain the corresponding T- a - 
Pareto optimal solution, say x*. 

5) To test whether this T- a -Pareto optimal solution x* to (3) is also a -Pareto optimal solution 
to (2), solve (4). 

6) Let x and £ be optimal solution of (4) in step 5. Then two cases may arise: 
a) If e i = 0 V i , then x* is a -Pareto optimal solution to a -MOLPP. 
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b) Iff. > 0 for at least one j , then x* is not <2 -Pareto optimal solution to <2 -MOLPP. Instead 
of x*, x is a -Pareto optimal solution to the a -MOLPP. 

7) This solution is T- a -Pareto optimal solution as well as a -Pareto optimal solution for fixed 
a(0 < a < 1) to the a -MOLPP. And the algorithm is complete. 

VI. Numerical Examples 

A. Numerical Example 1 . 

Here we consider an imprecise multiobjective linear programming problem with fuzzy 

coefficients as follows: 

imprecise max z : = 0.875xj - C n x 2 
imprecise min z 2 = C 2l x { + x 2 
imprecise max z 3 = C 31 x : - 2x 2 
subject to 

5x, + 7x 2 < 12, 9x 1 +x 2 < 10, -5jc, + 3x 2 < 3, x l + x 2 > 1, 
x { > 0, x 2 > 0. 

For the sake of simplicity alone, we denote 

X ={(x l ,x 2 ):5x l +7x 2 <12, 9x t +x 2 < 10, -5x;+3x 2 <3,x l +x 2 > 1, x l > 0, x 2 > 0} 

Here the coefficients C 12 , C 21 , C 31 are fuzzy numbers with membership functions given by 

jU Cn (C 12 ) = max(l-|C 12 -l|,0), 
jU Cii (C 21 ) = max (l - |C 21 - 5| , 0) , 
M Cn (C 31 ) = max (l - 2 |C 12 - 6| , 0) . 

Suppose that the DM choses value of the parameter as a = 0.7. Then for a = 0.7 , the fuzzy 
coefficients C n ,C 2l ,C 3l are triangular fuzzy numbers as 
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C 12 =[0.7, 1.3], C 21 =[4.7, 5.3], C 31 =[5.85, 6.15]. We first solve this problem by using classical 
membership function of fuzzy set theory. 

1 ) Imprecise a -MOLPP with first set of parametric T -characteristic functions analogous to 
usual membership functions: 

To construct membership functions for each of the imprecise objectives, we find individual 
maximum, and minimum of each of the objective functions (properly using 0.7-cut of triangular 
fuzzy coefficients) under the given constraints as follows: 

TABLE I 

Individual Maximum, and Minimum 
Objective Functions Individual Maximum Individual Minimum 

Zl U) = 0.875*! -1.3* 2 0.97 - 1.7 

z 2 (x) = 4.7x 1 +x 2 5.7 1 

z 3 (jt) = 6.15*! -2;t 2 6.8 -2 



Next suppose that the DM specifies goal, and tolerance for each of the imprecise objective 
functions as follows: 

TABLE II 

Goals, and Tolerances of Objective Functions 



Objective Functions 



Goals 



Goal plus Tolerance 
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Values 

zi(x) 09 - 1 

z 2 (x) 1.5 3 

z 3 (x) 6.5 3 



Based on the assigned goals, and tolerances, the membership functions of each of the 
imprecise objective functions may be formulated as follows: 

[ 1 , if 0.875^ -1.3x2 > 0.9 

/ v i \ 0.875x -1.3x, +1 
W) =^M = \ g ^ . </ -1^-875^ -1.3x 2 <0.9 

[ 0 , if 0.875^ -1.3x2 <-l 



r 2 (c 2 x)^(c 2 x) = 



, ;/4.7x 1 +x 2 <1.5 



3-(4.7x 1 +x 2 ) 



if 1.5<4.7x 1 +x 2 <3 
( /4.7x 1+ x 2 >3 



r 3 (c 3 x)=^(c 3 x): 



16.15^-2x2-3 
: 15 ' 



,i/6.15x 1 -2x 2 >6.5 

if 3<6A5x l -2x 2 <6.5 
if 6.15^-2*2 < 3 



If the DM does not specify goals, and tolerances, individual maximum, and minimum of 
objectives may be used to construct membership functions. 

Then imprecise MOLPP with fuzzy coefficients becomes a single objective linear 
programming problem as follows: 
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max A 
subject to 
0.875x, - 1.3x 9 +1 




0 < A < l,x l ,x 2 > 0. 



Using Lingo 15.0.32, we find that this problem has no feasible solution. 

2) Imprecise a -MOLPP with second set of parametric T -characteristic function analogous to 
redefined membership functions: 

Next we construct redefined membership functions //..(c,.*), 2 = 1,2,3 (as proposed by Wu et al 
[17]) based on the goal, and tolerance (assigned by DM) of each of the imprecise objectives as 
follows: 



0.875^-1.3x2+1 



if 0.875x 1 -1.3x 2 >-l 
if 0.875x 1 -1.3x 2 <-l 



1.9 
0 



T 2 (c 2 x)= J U R 2 (c 2X ) 



' 3-(4.7x 1+ x 2 ) ^ 



-, /f4.7x 1 +x 2 <3 
, if 4.7x 1 +x 2 >3 



1.5 
0 




Based on these redefined membership function, we convert imprecise MOLPP with fuzzy 
coefficients to a single objective linear programming problem [4] as follows: 
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max X 
subject to 
0.875^-1.3*2+1 
L9 

3 -K + * 2 L , 

1.5 

6.15x -2x 0 -3 



5^ +lx 2 < 12, 9x { +x 2 < 10, -5^+3*2 < 3,x { + x 2 > 1, 
X,Xy,x 2 >0. 

Using Lingo 15.0.32 to solve this problem, we find that the problem has no feasible solution. 

3) Imprecise a -MOLPP with third set of parametric T- characteristic functions: 
Next we introduce another set of parametric T-characteristic functions for each of the 

imprecise objective functions such that each of the parametric T-characteristic functions are 

strictly monotonic over entire domain and higher membership grade of parametric T- 

characteristic functions give better result for objectives. 
Based on goals, and tolerances of imprecise objective functions, we define parametric T- 

characteristic functions for each of the imprecise objective functions as 

0.875x 1 -1.3x 2 +l 3-(4.7x.+x 2 ) 6.15X.-2X.-3 
T A c i x )= Y9 ,V(x„x 2 )eX, T 2 (c 2 x)= — ,V(x 1 ,x 2 )eX, r 3 (c 3 x)= — ,V(x 1 ,x 2 )eX. 

Here membership grade between zero , and one represents partly presence. Membership grade 
more than one represents must presence whereas membership grade less than zero represents 
never presence. Using these parametric T-characteristic functions, we convert imprecise MOLPP 
with fuzzy coefficients to single objective linear programming problem as follows: 
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max X 
subject to 

0.875^-1.3*2+1 3-(4.7x 1 + x 2 ) 6.15^-2*2-3 

> X, > X, > X, 

1.9 1.5 3.5 

5xj + 7*2 < 12, 9x^ + *2 < 10, - 5x^ + 3x2 < 3, x^ + X2 > 1, Xp X2 > 0, X is unrestricted in sign. 

Using Lingo 15.0.32, we obtain T- a- Pareto optimal solution as 

A* = -0.0873883, x,* = 0.5759682, x 2 * = 0.42403 18, Zl (x*) = -0.047269 12, z 2 (x*) = 3.131082, 
z 3 (x*) = 2.694141, T x (z, (x*)) = 0.5014373, r 2 (z 2 (x*)) = -0.08738828, T 3 (z 3 (x*)) = -0.08738828. 

Then we test a- Pareto optimality of T- a- Pareto optimal solution by solving the following: 



subject to 

(0.875x, -1.3x 2 ) + l 

^ ^ e, > 0.5014373, 

1.9 

3~(4.7x 1+ x 2 ) 



e 2 > -0.08738828, 
£>-0M 



1.5 

(6.15^-2x^-3 
15 

5xj+7x 2 <12, 9x 1 +x 2 <10, -5xj+3x 2 <3,Xj+x 2 >1, 
x p x 2 , f 2 , e 3 >0. 

Using Lingo 15.0.32, we obtain a- Pareto optimal solution x and £ as follows: 

e x =e 2 =e 3 =0, x, =0.5759682, x 2 = 0.4240318, z,(x) = -0.04726912, z 2 (x) = 3.131082, z 3 (x) = 2.694141, 
7; ( Zl (x)) = 0.5014373, T 2 (z 2 (x)) = -0.08738828, T 3 (z 3 (x)) = -0.08738828. 

Here all e t = 0\/i = 1,2, 3. Hence the solution obtained in previous stage is itself a -Pareto 
optimal solution to a-MOLPP. 



4) Imprecise a -MOLPP with fourth set of parametric T- characteristic functions: 
Again we may construct parametric T-characteristic functions for each of the imprecise 
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objective functions as follows: 

T 1 (c i x) = 0.i75x l -l.3x 2 Mx l ,x 2 )e X, T 2 (c 2 x)=-{4Jx 1 +x 2 ),\/(x 1 ,x 2 )e X, 7 3 (c 3 x) = 6.15*, -2x 2 ,V(x 1 ,x 2 )e X. 
Here the parametric T-characteristic functions are again strictly monotonic over entire domain, 

and higher membership grade of parametric T-characteristic functions give better result for 

objectives. But membership grade zero or one of parametric T-characteristic functions do not 

represent grade of membership in usual fuzzy sense. Precisely we have picked up the only 

necessary property (strict monotonicity over entire domain) of membership function to construct 

our proposed parametric T-characteristic functions. 

Based on these parametric T-characteristic functions, we convert the a-MOLPP into single 

objective linear programming problem [5] as follows: 

max X 
subject to 

0.875xj —1.3;t 2 > X,—{4.7x l + x 2 )> X, 
6A5x 1 -2x 2 > X, 5x l +lx 2 < 12, 
9x l +x 2 <\0,-5x l +?>x 2 <3, 
x x +x 2 > 1,.*,, x 2 > 0,X is unrestricted in sign. 

Using Lingo 15.0.32, we obtain the T- a- Pareto optimal solution as 

A* = -1.312236, *j* = 0.08438819, x 2 * = 0.91561 18, Zl (x*) = -1.116456, z 2 (x*) = 1.312236, z,(x*) = -1.312236, 
7^ (zj (jc*)) = -1.1 16456, T 2 (z 2 (x*)) =-1.312236, T 3 (z,(x*)) = -1.312236. 

To test a -Pareto optimality of this T- a -Pareto optimal solution, we solve the following 
problem: 
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3 

max ^ e i 

(0.875*, -1.3* 2 )-£, > -1.116456, 

-(4.7*! + x 2 )-e 2 > -1.3122363, 

(6.15*, -2* 2 )-£ 3 > -1.3122363, 

5*, + 7x 2 < 12,9*, + * 2 < 10,-5*, + 3* 2 < 3,*, + * 2 > 1, 

*,, * 2 , e 2 , e 3 > 0. 

Using Lingo 15.0.32, we obtain the a- Pareto optimal solution * and e as 

e l= e 2 =e 3 =0, I = -1.312236, 

3t| =0.08438819, x 2 =0.9156118, z,( x) = -1.116456, z 2 (x) = 1.312236, z 3 (x) = -1.312236, 
T x ( Z[ (j)) = -l.l 16456, T 2 (z 2 (x)) = -1 .3 12236, T 3 (z 3 (*)) = -1 .3 12236. 

Here all =0V/ = 1,2,3. Hence the solution obtained in first stage x* is itself or -Pareto optimal 
solution to the a- MOLPP. 

Therefore usage of membership function or redefined membership function gives no feasible 
solution to the problem. But if we use well defined parametric T-characteristic functions, we may 
find a -Pareto optimal solutions. If we compare results obtained by using first set of proposed 
parametric T-characteristic functions with usual membership functions, we observe that the 
membership grade of zi(x) is more than 0.5; certainly much more acceptable scenario than 
infeasible solution to some problem. Moreover second set of parametric T-characteristic 
functions gives another a- -Pareto optimal solution where conventional membership grade of 
objective function Z2(x) is 1, and redefined membership grade is 1.12. 

TABLE III 

Status of Solutions in Different Cases 



Parametric T-characteristic Function 



Status 
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of Solution 


First set: Analogous to usual membership functions 


No feasible solution 




Second set: Analogous to redefined membership 
functions 


No feasible solution 




Third set: Usage of strictly monotonic part of 


a-Pareto optimal 


solution 


membership function over entire domain 


exists 




Fourth set: Usage of any strictly monotonic linear 


a-Pareto optimal 


solution 


function 


exists 





B. Numerical Example 2: 



We consider another multiobjective linear programming problem on investments made by 
private bank, based on paper by J. L. Eatman, and C. W. Sealey Jr. [8]. To match with reality, we 
may assume that DM has imprecise goal for each of three objectives: a) maximization of profit, 
b) minimization of capital adequacy ratio, and c) minimization of illiquid risk assets. Also 
suppose that some coefficients viz. C 15 ,C I6 ,C 25 ,C 28 ,C 31 are triangular fuzzy numbers. So the 
MOLPP now becomes MOLPP with fuzzy coefficients under imprecise environment. 

Then linear programming problem may be formulated as follows: 
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imprecise max 0.04x 2 +0.045x 3 +0.055x 4 +C l5 x 5 +C l6 x 6 +0.085x 7 +0.092x g 
imprecise min ^ (0.005x 2 + 0.04x 3 + 0.05x 4 + C 25 x 5 +0. lx 6 + 0. lx 7 + C 28 x 8 ) 

imprecise min -^(C 3l x 6 + x 1 + x g ) 
subject to 

x, + x 2 +x 3 + x 4 + x 5 + x 6 + x 7 +x & =(20+150+80), 
jc 1 >(0.14x150)+(0.04x80), 

x, +0.995x 2 +0.96*3 +°- 9 *4 +0.85x 5 > (0.47xl50) + (0.36x80), 
Xj > 0.05(20+150+80), Vj =1...8, 
jc 8 > 0.3(20+150+80),^...^ >0 

Here C 15 ,C 16 ,C 25 ,C 28 ,C 31 are all fuzzy coefficients in the form of triangular fuzzy numbers as 
C 15 =(0.06A07A08),C 16 =(0.10A105A110),C 25 =(0.070A075A080),C 28 =(0.095,0.10,0.105),C 31 =(0.8,1.0,L2). 
We suppose that the DM supplies aspiration level a = 0.7 for these fuzzy coefficients. 

1 ) First set of parametric T -characteristic functions: 

Here we take T-characteristic functions for each of imprecise objective functions analogous to 
membership functions of fuzzy set theory. And each of the parametric T-characteristic functions 
is strictly monotonic over entire domain and higher membership grade of parametric T- 
characteristic functions give better result for objectives. Accordingly, we need to determine the 
goals, and tolerances of each of the imprecise objective functions. Therefore we compute 
individual maximum, and minimum of each of the imprecise objective functions under the given 
constraints and obtain the results as follows: 

TABLE IV 
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INDIVIDUAL MAXIMUM, AND MINIMUM 



Objective Function 



Individual Individual 
Maximum Minimum 



z l (x) = c ] x 
z 2 (x)=c 2 >c 
z 3 (x)=c 3 x 



18.89438 11.95625 
0.9338129 0.5996875 
7.313125 4.9625 



Next, suppose that the DM specifies the goal, and tolerance for objective function cjx only. For 
other objective functions, we use individual maximum, and minimum values of objectives from 
table no 4 and construct parametric T-characteristic functions T t (c ; x), i = 1,2,3, analogous to 
membership functions of fuzzy set theory as follows: 



x-U 0.9338129-c,x 7.313125-c-j 
TAc,x) = ^ VxeX, r o (c 0 x) = ^-VieX, T~(c~x)= ^-Viel. 



0.3341254 3 J 2.350625 

Based on these parametric T-characteristic functions, we convert the a-MOLPP into single 

objective linear programming problem [5] as follows: 

max 1 
subject to 

-(0.04x 2 +0.045x,+0.055x.+0.07 3x 5 +0.1065x 6 
2 

+ 0.085x 7 +0.092x 8 - 14) > X, 

{0.9338 129 - 0.05(0.005 x 2 + 0.04x 3 + 0.05x 4 

0.334 1254 2 3 4 

+ 0.07 35 x 5 + 0.1x 6 + 0.1x 7 + 0.0985 x 8 )} > X, 



1 



(7.3 1 3 125 - 0.05 (0.94x 6 + x 7 + x 8 )} > X, 
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x l + x 2 + x 3 + x 4 + x 5 +x 6 +x 1 + x 8 =(20 + 150 + 80), 
x, >(0.14xl50)+(0.04x80),x. > 0.05(20+150+80), forj =1...8, 
x, + 0.995x 2 +0.96x 3 +0.9x 4 +0.85x 5 >(0.47xl50) + (0.36x80), 
x 8 > 0.3(20 + 150 + 80), A is unrestricted in sign, x l . . . x 8 > 0. 

Using Lingo 15.0.32, we obtain the T- or-Pareto optimal solution as 

A* = 0.81,x 1 * = 24.2, x 2 * = 78.84,x 3 * = x 4 * = x 5 * = x 7 * = 12.5,x 6 * = 21.96, x s * = 75, Zl (x*) = 15.62, z 2 (x*)=0.66, 
z 3 (x*) = 5.41, ( Zl (x*)) = T 2 (z 2 (x*)) = 0.8 1, T 3 (z 3 (x*)) = 0.8 1 1 . 

To test a-Pareto optimality of this T-or-Pareto optimal solution, we solve the following 
problem 

max ^ £ i 
subject to 

y(0.04x 2 +0.045x 3 +0.055x 4 +0.073x 5 +0.1065x 6 +0.085x 7 +0.092x 8 -\A)-e l > 0.8087605, 

{(0.9338 129 -0. 05(0. 005x 2 + 0.04x 3 

0.3341254 2 

+0.05x 4 + 0.0735x 5 + 0.1x 6 + 0.1x 7 + 0.0985x 8 )} - £ 2 > 0.8087605, 

{7.3 13 125 - 0.05 (0. 94x 6 + x 7 + x 8 )} - e, > 0.8107958, 

2.350625 v 6 7 8 /13 

X; + x 2 + x 3 + x 4 + x 5 + x 6 + x 7 + x 8 = (20 + 150 + 80), Xj > (0. 14 x 1 50)+(0.04 x 80), 

x x + 0.995x 2 + 0.96x 3 + 0.9x 4 + 0.85x 5 > (0.47 x 150) + (0.36 x 80), 

x. > 0.05(20+150+80), for j=l,2...8,x 8 > 0.3 (20 + 150 + 80), Xl ...x s ,e x ,e 2 ,e } > 0. 

Using Lingo 15.0.32 to solve this problem, the a-Pareto optimal solutions x and fare as 
follows: 

e x = s 2 = e 3 = 0, X = 0.81, jc, = 24.2, x 2 = 78.84, x 3 = x 4 = x 5 = x 7 = 12.5, x 6 = 21.96, % = 75, 
z, (x) = 15.62, z 2 (x) = 0.66, z 3 (x) = 5.41, T x (z, (x)) = T 2 (z 2 (x)) = 0.8 1 T 3 (z 3 (x)) = 0.8 1 1 . 

Here all £ t = 0V/ = 1,2, 3. Hence the solution obtained in first stage i.e. x* is itself a-Pareto 
optimal solution to thea-MOLPP. 
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2) Second set of parametric T -characteristic functions: 

Instead of using parametric T-characteristic function analogous to membership function of 
fuzzy set theory, we construct another set of parametric T-characteristic functions to obtain a- 
Pareto optimal solution of or-MOLPP for fixed a =0.7 as follows (so that each of the parametric 
T-characteristic functions are strictly monotonic over entire domain , and higher membership 
grade of parametric T-characteristic functions give better result for objectives) 

7] (c,jc) = 0.04x 2 + 0.045x 3 + 0.055x 4 + 0.073x 5 + 0.1065x 6 + 0.085x 7 + 0.092x 8 , Vxe X 
T 2 (c 2 x) = -^(0.005x 2 +0.04x 3 +0.05x 4 +0.0735x 5 +0.1x 6 +0.1x 7 + 0.0985x 8 , Vxe X 
r 3 (c 3 x) = -^(0.94x 6 +x 7 +x 8 ), Vxe X 

Based on these parametric T-characteristic functions, we now convert the or-MOLPP into 

single objective linear programming problem [6] as follows 

max X 
subject to 

0.04x 2 +0.045x 3 +0.055x 4 +0.073x 5 +0.1 065x 6 +0.085x 7 +0.092x 8 > A, 

-0. 0 5 (0. 005 x 2 + 0.04 x 3 + 0.05 x 4 + 0.07 3 5 x 5 + 0.1 x 6 + 0.1 x 7 + 0.09 8 5 x 8 ) > X, 

-0.05 (0.94x 6 + x 7 + x 8 ) > X, 

X! + x 2 + x 3 + x 4 + x 5 + x 6 +x 7 + x 8 = (20 + 150 + 80), 
x x > (0.14 x 150)+(0. 04 x 80), 

x l + 0.9 9 5x 2 + 0.9 6x 3 + 0.9 x 4 + 0.8 5 x 5 > (0 .47 x 1 5 0 ) + (0 .3 6 x 8 0 ) , 

x. > 0.05(20+ 150+80), forj=1...8, 

x 8 > 0.3 (20 + 1 50 + 80), 

Xj ... x 8 > 0,X is unrestricted in sign. 

Using Lingo 15.0.32, we obtain the T- or-Pareto optimal solution as 

X* = -4.96, x^ = 24.2, x 2 * = 88.3, x 3 * = x 4 * = x 5 * = x 6 * = x 7 * = 12.5, x 8 * = 75, 

Zi(x*) = 14.99, z 2 (x*) = 0.62, z 3 (x*) = 4.96, T { (z t (x*)) = 14.99, 7/ 2 (z 2 (x*)) = -0.62, T 3 (z 3 (x*)) = -4.96. 
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To test a-Pareto optimality of this T-or-Pareto optimal solution, we solve the following 
problem 

3 

max £ £■ 
i=\ 1 

subject to 

0.04x 2 +0.045x 3 +0.055x 4 +0.073x 5 +0.1065x 6 +0.085x ? +0.092x g - ^ > 14.98825, 

-0.05 (0.005x 2 +0.04x 3 +0.05jc 4 +0.0735x 5 +0.1x 6 +0.1x ? + 0.0985x g - e 2 > -0.6186375, 

-0.05 (0.94x 6 + Xj + x g ) - e 3 > -4.9625, 

x l + x 2 + jc 3 + x 4 + x 5 + x 6 + x 7 + x g = (20 + 150 + 80), 

Xj > (0.14xl50)+(0.04x80), 

Xj +0.995x 2 + 0.96x 3 +0.9x 4 +0.85x 5 > (0.47 x 150) + (0.36 x 80) , 

xj > 0.05(20+150+80), for j =1...8, 

x 8 > 0.3(20 + 150 + 80), x v ..x 8 ,£ v e 2 ,£ 3 > 0. 



Using Lingo 15.0.32 to solve this problem, we obtain the a-Pareto optimal solution x and £ 



e x = e 2 = £ 3 = 0, X = -4.96, ^ = 24.2, x" 2 = 88.3, x" 3 = x" 4 = x" 5 = x^ = x" ? = 12.5, x g = 75, 

Z j(x) = 14.99, z 2 (x) = 0.62, ^(x) = 4.96, 2^ (3e)) = 14.99, T 2 (z 2 (x)) = -0.62, = -4.96. 

Here all e i =0V/ = l,2,3.Hence the solution obtained in first stage x* is itself a-Pareto 
optimal solution to the a- MOLPP for fixed a . 



TABLE V 

Different Pareto Optimal Solutions 



Optimal Values Using 



Objective Functions 



First Set of Parametric Second Set of Parametric 

T-characteristic Functions T-characteristic Functions 
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Profit 



15.62 



14.99 



Capital-adequacy Ratio 



0.66 



0.62 



Illiquid Risk Assets 



5.41 



4.96 



Clearly first set of parametric T-characteristic functions in a-MOLPP results in higher profit 
whereas second set of parametric T-characteristic functions gives better result for capital- 
adequacy ratio as well as illiquid risk assets, two important pillars for optimal decision of this 
example. And both are or-Pareto optimal solutions. 



We have observed that Wu et al (2015) [17] did not use redefined membership function 
according to their proposed definition. In numerical example one of this paper, we have shown 
that although redefined membership function in ar-MOLPP correctly identifies the best a- 
Pareto optimal solution among the good cases; it fails to identify bad a-Pareto optimal solution 
among the worst cases! And proper usage of membership or even redefined membership 
functions may suggest an a-MOLPP infeasible, although we may find or-Pareto optimal 
solution(s) by using parametric T-characteristic functions. 

On the other hand, in real life scenarios, DM may not be able to identify the goal or tolerance 
or both; and in some cases, DM may not be present at all. In these cases, it is difficult to 
construct membership function or redefined membership function in imprecise environment. 
And the usage of well-defined parametric T-characteristic function decreases frequent interaction 
with DM [11], and thus it saves precious time of both. 



VII. CONCLUSIONS 
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Further higher value of membership function generating better value for objective is one of the 
primary aims of using membership function in imprecise environment. And if we use usual 
membership function or redefined membership function to solve MOLPP with fuzzy 
coefficients, we can find at most one a-Pareto optimal solution in imprecise environment. 

But if we use parametric T-characteristic functions, we may find different solutions that are all 
cr-Pareto optimal. Finally DM is the king, and his choice is always final. And scientific, and 
especially optimization tools must try to supply more than one 'best choice' to DM, if possible; 
and that should hold even more in imprecise environment. So usage of parametric T- 
characteristic function in solving MOLPP with fuzzy coefficients is more suitable to DM than 
membership or redefined membership function under fuzzy environment. 

Hence in a -MOLPP under imprecise environment, we propose to take general procedure of 
using parametric T-characteristic functions covering the entire set of possible solutions in place 
of membership functions and redefined membership functions. From all these we may add that 
under these circumstances the issue of infeasibility in a -MOLPP even by using fuzzy set theory 
yielded from several published methods seems worthwhile or even necessary to reconsider. 
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